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I - THERMODYNAMIOUE DES MACHINES DE COMPRESSION ET DE DETENTE : 



Les machines de compression et de detente peuvent etre schematisees suivant les deux grands 
modeles de la figure 1 oil l'on peut noter les grandeurs les plus significatives : 

• Le debit masse q et les vitesses absolues Ci et C2. 

• Les caracteristiques d'etat du fluide a l'entree (pi, Ti) et a la sortie (p2, T 2 ). 

• La chaleur regue par le fluide Qe entre l'entree et la sortie. Par convention, Qe est comptee 
positivement si elle est regue par le fluide et negativement dans le cas contraire. 

• T t h le travail echange entre le fluide et l'operateur de la machine (piston, palette, roue 
ailetee, etc....). On l'appelle egalement : travail indique, travail moteur, travail interne,... 

Pour une machine de compression x t h est regu par le fluide, pour une machine de detente, c'est 
le fluide qui communique x t h a l'operateur. 

Le travail de frottement Xf est directement lie a la viscosite du fluide. Pour les fluides 
compressibles de faible masse volumique, le travail des forces de pesanteur est generalement 
negligeable. 




Compression Detente 

Figure 1 : Machines de compression et de detente. 



Suivant le cas, le travail sur l'arbre est superieur ou inferieur au travail theorique. Les pertes 
mecaniques x m materialisent les frottements sur les paliers et le dispositif d'etancheite. 
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1-1 - Expressions du travail indique : 



Le premier principe de la thermodynamique, ecrit pour un systeme ouvert en mouvement 
permanent, est utilise pour les machines a fluide : 




( 1 ) 



avec A : valeur finale - valeur initiale ou valeur de sortie - valeur d'entree. 
+ : machine de compression. 

- : machine de detente. 



L' enthalpie h du fluide est introduite par : 



h = J^ + jTdS (2) 

C 2 

l'enthalpie totale H : H = h + — — (3) 

l'entropie S : TdS = d (Qe+Xf) (4) 

Ainsi, on retiendra deux formes principals du travail indique : la forme energetique (5) 
issue directement de (1) et la forme dynamique (6) issue de (1), (2) et (4) : 



T th =± 


C 2 

Ah + A Qe 

2 




rl 

s- 

II 

1+ 


r 2d P A c 2 
b A b Tf 

p 2 1 





(5) 

(6) 



avec + : machine de compression 
- : machine de detente 

L'equation dynamique (6) ne peut etre integree que dans quelques cas particuliers ou 
1'evolution de la masse volumique du gaz est connue et facilement modelisable comme nous 
le verrons ulterieurement. 

Une derniere forme est egalement utilisable pour les turbomachines (machines 
rotodynamiques ou les echanges energetiques s'effectuent principalement par variation de la 
quantite de mouvement du fluide entre l'entree et la sortie du rotor), c'est la forme cinematique 
(ou equation d'Euler des turbomachines) deja examinee : 



T th = ±[U A Cu\ 



(7) 



± memes conventions de signe. 
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Equation d'etat des gaz parfaits 



Elle permet de calculer la troisieme variable d'etat a partir des deux autres, nous retiendrons la 
loi des gaz parfaits sous sa forme etablie par Mariotte : 



avec r : constante du gaz. 

On demontre que pour les gaz parfaits, on peut ecrire 

dh = Cp dT 

Suivant la nature du gaz, Cp et Cv varient faiblement ou non en fonction de la temperature 



(tableau 1). 






Runnels : 


Cp 
Y = — 
Cv 


(10) 




r = Cp - Cv 


(11) 




Cp= 

Y-l 


(12) 




Cv = r 
Y-l 


(13) 



( 8 ) 



(9) 



Tableau - Valours, pour differents gaz, de la masse molaire, de la constante r 
et de c p (en J ■ kg -1 • K -1 ) en fonction de la temperature 


I 

1 

Nature du gaz 


Symbole 


Masse 

molaire 

(en 

kg • mol' 1 ) 


r 

(en 

J ■ kg-’K-'l 


0 


15 


100 


200 


300 


Temp 

(e 

400 


6rature 

n“C) 

500 


600 


700 


800 


900 


1 

1000 


Air 

Azote 

Ammoniac 

Dioxyde de carbone .... 

Dioxyde de soufre 

Ethylene 

Hydrogfcne 

Methane 

Monoxyde d'azote 

Monoxyde de carbone 

Oxyde d'azote 

Oxyg6ne .. 

Sulfure d'hydrogfcne .. 

Vapeur d'eau loin 
de la liquefaction 


n 2 

nh 3 

co 2 

so 2 

c 2 h 4 

h 2 

ch 4 

NO 

CO 

n 2 o 

0 2 

h 2 s 

h 2 o 


0,02898 

0,028 

0,017 

0,044 

0.064 

0,028 

0,002 

0,016 

0,030 

0.028 

0,044 

0,032 

0,034 

0,018 


287 

297 

489 

189 

130 

297 

4157 

520 

277 

297 

189 

260 

244 

462 


1004 

.• 

1037 

2054 

820 

606 

1410 

14229 

2155 

997 

1039 

849 

916 

1025 

1853 


1005 

1037 

2078 

835 

616 

1516 

14279 

2190 

997 

1041 

865 

916 

1030 

1856 


1014 

1046 

2205 

920 

665 

1854 

14438 

2452 

997 

1046 

954 

933 

1062 

1881 


1025 

1052 

2387 

996 

715 

2201 

14521 

2795 

1012 

1058 

1025 

962 

1111 

1933 


1040 

1069 

2573 

1067 

753 

2498 

14542 

3172 

1033 

1079 

1083 

996 

1163 

1987 


1070 

1092 

2736 

1121 

786 

2749 

14634 

3498 

1062 

1104 

1129 

1025 

1213 

2054 


1092 

1117 

3931 

1163 

807 

2971 

14660 

3812 

1088 

1129 

1175 

1050 

1268 

2117 


1113 

1138 

3096 

1201 

823 

3151 

14773 

4084 

1109 

1159 

1209 

1075 

1318 

2184 


1134 

1163 

3255 

1230 

841 

3327 

14940 

4321 

1129 

1180 

1242 

1080 

1364 

2255 


1155 

1182 

3398 

1255 

849 

3477 

15128 

4540 

1146 

1201 

1263 

1100 

1406 

2326 


1169 

1201 

3532 

1270 

857 

3607 

15317 

4758 

1159 

1217 

1280 

1113 

1439 

2393 


1184 

1215 

3653 

1276 

866 

3724 

15526 

4942 

1171 

1237 

1297 

1123 

1473 

2456 



Pour fair, on retiendra generalement : r = 


287 


J/kg°K 


Y 


1,4 




Cp = 


1000 


J/kg°K 


Cv = 


713 


J/kg°K 
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Cas particuliers : 



Quelques cas particuliers permettent une application rapide et tres explicite de ces equations 
generates : chaudiere, echangeur, chambre de combustion. Dans ces appareils, pas 
d'echange de travail moteur avec l'exterieur et variation d'energie cinetique negligeable ; la 
forme energetique donne : 



Ah = Qe = J 2 CpdT 



(14) 



soit pour Cp constant : Ah = Cp (Ti - Ti) 

laminage dans une vanne calorifugee : 



Xth = Qe = 0 



Lorsque la variation d'energie cinetique n'est pas prise en compte, l'equation (5) se reduit a : 



Ah = 0 ou h = cte 



(15) 



soit une evolution isenthalpique. 
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1-2 - Compression et detente isentropiques : 



C'est 1'evolution retenue pour les fluides compressibles dans le cas particulier reunissant 
simultanement les deux conditions : 



• evolution adiabatique : Qe = 0 

• fluide non visqueux : Xf = 0 

Soit encore d'apres (4) : 



T dS = 0 ou S = cte 



L'equation de Laplace est l'equation retenue pour decrire cette evolution isentropique : 




(16) 



y : exposant isentropique defini en (10). 



Evolution de la temperature : 



En retenant 1 comme indice initial et 2 comme indice final, 1'evolution de la temperature est 
calculable par une combinaison de (16) et (8). 





y-i 

y 



(Pi) 



(17) 



L'indice 2S caracterise 1'evolution isentropique (figure 2). Les temperatures sont exprimees en 
°K et les pressions sont des pressions absolues. 




Figure 2 : Compression et detente isentropiques dans le diagramme T-S. 
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Travail moteur : 



L'expression du travail moteur est calculable en exprimant p a partir de l'equation de Laplace 
(16). 

Pour 1'evolution isentropique, on obtient d'apres (6) : 

7-1 
7 

et d'apres (17) et (12) : 




On pouvait obtenir directement cette equation (19) a partir de la forme energetique du travail 
moteur (5) et de l'equation des fluides parfaits (9) pour Cp constant. 

( 20 ) 

Enfin, en negligeant la variation d'energie cinetique, le diagramme h-S donne directement le 
travail moteur. Pour Cp constant, les isothermes sont des horizontales (figure 3). 







Compression 



Detente 



Figure 3 : Compression et detente isentropiques dans le diagramme h-S. 



MACHINES A FLUIDES COMPRESSIBLES - TURBINES ET COMPRESSEURS 

Robert REY - Professeur 



7 






1-3 - Evolution polytropique : 



En fluide reel, les frottements internes if se dissipent en chaleur et conduisent, en fin de 
compression ou de detente, a une temperature differente de celle correspondant a 1'evolution 
isentropique. A ce phenomene se superpose, en s'additionnant ou se retranchant, l'echange de 
chaleur Qe avec l'exterieur. La connaissance precise de 1'evolution reelle necessiterait une 
mesure continue de la pression et de la temperature en un grand nombre de points a l'interieur 
de la machine. Seules les donnees d'entree et de sortie etant generalement disponibles, on cale 
1'evolution sur une caracteristique de forme similaire a celle de l'isentropique. On definit ainsi 
un exposant polytropique k a partir de la relation : 



h 

T, 



f \ 

Pi 

P iJ 



k-1 



( 21 ) 




Compression Detente 

Figure 4 : Evolution polytropique dans le diagramme T-S. 

Pour une machine calorifugee (evolution adiabatique) les mesures effectuees a l'entree et a la 
sortie de la machine permettent le calcul de k. 

k>l,4 pour la compression adiabatique ^2) 

k<l,4 pour la detente adiabatique 



Rendement indique : 



C'est le rapport des travaux a foumir ou rcyu de la part du fluide par rapport a celui 
correspondant a 1'evolution isentropique : 

On aura d'apres (6) : 



*ls = 





r 


- dn C 2 




+ 




— + A + x r 




J 


P 2 




± 


[ 2 — + aC-" 






Jl p 2 





(23) 
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En negligeant l'energie cinetique, le rendement indique se transforme en : 



Pour la compression : 



Pour la detente : 



Cp(T ;s .T,) 
% Cp(T 2 .T,) 



(P 2 /P 1 ) l 

k-1 

(P 2 /Pl)~ - 1 



_ Cp (T 2 T, ) (pp/pj'k'-l 
-Cp(T 2S .T 1 )- (p2/pi) ^- 1 



( 24 ) 



(25) 



Ces grandeurs sont calculables a partir des mesures a l'entree et a la sortie. 

Lorsque P 2 /P 1 est proche de 1, un developpement limite au premier ordre, donne les valeurs 
approchees : 



Compression : 



Tls 



y-1 k 



y k-1 



(26) 



Detente 



^8 = 



k-1 y 
k y-1 



(27) 



Ces expressions permettent de lier directement l'exposant polytropique k avec le rendement 
indique. Elies sont particulierement interessantes au stade de l'avant projet. Dans le 
diagramme h-S le rendement indique est directement accessible. 




Compression Detente 

Figure 5 : Evolution polytropique dans le diagramme h-S. 



Compression : 


Ahs 
s Ah 


(28) 


Detente : 


Ah 
s Ahs 


(29) 
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II - TURBOMACHINES DE DETENTE. 



Les applications des turbines a gaz et a vapeur sont nombreuses et variees. Elies sont 
justifiees par la grande puissance massique de ces machines : 

• Applications aeronautiques : turboreacteurs et turbomoteurs (Annexe 1 - Figures A1 a 
A4). 

• Applications spatiales : turbopompes. 

• Production electrique centralisee par unites de plusieurs centaines de MW (Annexe 1 - 
Figure A5). 

• Production electrique decentralisee avec ou sans cogeneration et/ou cycles combines 
(Annexe 1 - Figure A6). 

Fes enjeux economiques et industriels en relation avec ces machines sont done tres importants 
et relevent largement d'efforts et d'investissement multinationaux. 

Dans les turbines a gaz ou a vapeur, les etages sont de type a action ou a reaction. Fe degre de 
reaction caracterise la chute de d'enthalpie assuree par la roue ramenee a la chute d'enthalpie 
dans l'etage. Cette caracteristique est obtenue par une geometrie d'aubages tres caracteristique 
(Figure 6 et 7). 




0 1 2 



A - Etages a action 



f 


'A 


'E 


'// 






/z 


/z 























0 1 2 

B - Etages a reaction 



Figure 6 : Etages a action et a reaction [ref 1 ]. 
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Figure 7 : Formes typiques d'aubages rotoriques d action [1 ]. 

Pour conserver une vitesse debitante axiale constante, on observe une augmentation de 
l'envergure des aubages en cours de detente. Dans certains cas, l'ensemble est ceinture a la 
Peripherie (figure 8). En general, les aubages rotoriques sont rapportes sur le disque formant 
la roue (figure 9). 




Figure 8 : Turbine ci vapeur multistage a plan de joint horizontal 
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( Fives - Cail Babcock). 




assemblage pale/disque 
par pieds de sapin 

Figure 9 : Divers dispositifs de fixation des pedes [ref 2]. 



La majorite des turbines a vapeur sont a plan de joint horizontal avec des dispositifs 
d'etancheite sans contact de type labyrinthes (figures 10 et 11). 



Sfntnr 



Figure 10 : Dispositif d'etancheite vers 1'exterieur de type labyrinthe [ref 3]. 
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Les turbines a gaz sont constitutes de trois sous-ensembles distincts : les compresseurs, la 
chambre de combustion et les etages de turbines (figure 12). 

Suivant les cas, l’etage de detente est mono-arbre (turbine lite) ou separee en deux partie : 
l’etage ou les etages entrainant les compresseurs et les etages dits a turbine libre entrainant le 
recepteur ; altemateur, rotor d’helicoptere, etc... 

Les figures 12 et 13 montrent ces particularites : 




Entree d'air Compresseurs Chambre de Turbines Echappement 

combustion 



Figure 12 : Sous-ensembles constituant la turbine a gaz [2 ]. 
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R6ducteur 



Figure 13 : Turbomoteur ci turbine liee Makila 3G (Turbomeca). 

Le rendement global des turbines dependant directement de la temperature a la sortie de la 
chambre de combustion, la tenue mecanique des aubages du distributeur et du rotor du 
premier etage necessite un refroidissement. Ils sont aujourd’hui refroidis par un certain debit 
d’air preleve en sortie des compresseurs. Comme on peut l’observer sur la figure 14, la 
realisation des aubages en est rendue beaucoup plus difficile. 




Distributeur Roue 

Figure 14 : Ailettes refroidies. 
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2-1 - Configuration de l’etage de turbine : 



L’etage est constitue d’un distributeur (grille fixe) et d’un rotor tournant a la vitesse angulaire 
co. En phase de dimensionnement, on retient generalement les hypotheses suivantes : 

Le travail des forces de pesanteur est negligeable et la detente est adiabatique. 
le fluide est non visqueux et suit les lois des gaz parfaits. 

l’envergure des aubages est variable le long de l’axe d’ecoulement pour maintenir les 
vitesses debitantes constantes. 

Les vitesses absolues sont axiales a 1’ entree du distributeur et en sortie de roue. 
L’ecoulement est permanent et axisymetrique : l’ecoulement est moyenne dans les 
differentes sections d’ entree et de sortie des grilles fixes et mobiles. 



Une coupe cylindrique, effectuee a rayon r quelconque de l’ailetage et developpee en plan, 
donne Failure de la figure 15. L'hypothese d'un ecoulement bidimensionnel est un modele 
simplifie s'affranchissant des ecoulements secondaires et des couches limites (figures A7 et 
A8 de l'Annexe 1). La section 0 correspondant a l’entree du distributeur est a la pression p 0 , la 
section 1 correspondant a l’espace inter-grille est a la pression pi, la section 2 de sortie du 
rotor est a la pression p 2 . 



Espace statorique: 

Dans l’espace fluide statorique en fluide non visqueux (Xf = 0) l’equation de l’energie (6) 
s’ecrit pour x t h = 0 : 

H = f — + — = cte (40) 

J p 2 

ou encore pour la detente isentropique d’un gaz parfait, l’equation de Saint-Venant : 



H = Cp T + 



2 



Y P 

Y-l P 




(41) 




Espace statorique 



Espace rotorique 
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Figure 15 : Configuration de l'etage de detente 



Espace rotorique : 

Dans l'espace rotorique, en mouvement relatif, on introduit le terme de rothalpie qui, dans les 
memes conditions et hypotheses que precedemment, se maintient constante. 

En effet, en groupant les equations (6) et (7), on ecrira : 



1= [— + — — U Cu = cte 
J p 2 



(42) 



Par definition, la constante est notee I, c'est la rothalpie (reservee au rotor). D'apres les 
triangles de vitesses, elle s'ecrira egalement : 



I-f d P + w 2 


U 2 

— pfp 


J P 2 


2 



(43) 



Pour une machine axiale ou U = cte, l’analogie avec l’espace distributeur sera parfaite en 
remplagant la vitesse absolue par la vitesse relative. 



I = Cp T + 



w- 



y p 

1 ■ + = cte 



2 y-1 p 2 



L'equation de continuite s'ecrit ici par hypothese constructive : 

Ca 0 = Cai = Ca 2 - Ca 

Proprietes des triangles de vitesses : 

Ui = U 2 = U = cor 

Cui = Cu = Ca tg Pi + U = Ca tg ai 



(44) 



(45) 

(46) 

(47) 



U + Ca tg p 2 = 0 

Compte tenu des hypotheses faites, la variation d'enthalpie subie par le fluide a travers l'etage 
est calculable par l'equation d'Euler des turbomachines : 



Ahs = U . Cu 



(49) 



Degre de reaction. 

Le degre de reaction R est un terme de classification des etages de turbine, c'est le rapport 
entre la chute d'enthalpie dans le rotor et la chute d'enthalpie de l'etage complet figure 16. Soit 
en retenant les notations de la figure 15 : 
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( 50 ) 




Soit encore, d'apres les equations (42) et (45) et les hypotheses retenues : 

W 2 2 _ W, 2 

R = — ^ 

UCu 



En calculant W 2 etWi a partir des triangles de vitesses : 

W 2 2 = U 2 + Ca 2 2 
W! 2 = (Cu - U) 2 + Ca! 2 

on obtient finalement comme expression cinematique du degre de reaction : 




(51) 



(52) 




ACTION REACTION 



Figure 16 : Etage ci action et ci reaction. 
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Coefficients de pression et de debit. 

On introduit le coefficient adimensionnel de pression par la relation : 




( 53 ) 



La valeur de R reste au choix du constructeur. Une machine donnee sera capable de generer 
des chutes d'autant plus elevees que R sera petit. 

La limite inferieure de R est generalement conditionnee par la mise au point de profils 
satisfaisants (au sens des pertes) pour de tres grandes deflexions. 

La valeur de R = 0 est accessible dans de bonnes conditions avec des profils classiques 
(figure 7). 

On notera avec un grand interet que le degre de reaction nul se caracterise par une poussee 
axiale nulle sur la roue (pi = P 2 ). Par contre, dans les applications aeronautiques, la presence 
d'etages a reaction (par exemple R = 0,5) permet d'equilibrer la poussee axiale des 
compresseurs. 

Un second degre de liberte est necessaire pour fixer l'ensemble des triangles de vitesses. En 
general, on retient le coefficient de debit cp defini par le rapport : 




(54) 



La similitude des turbomachines (triangles de vitesses semblables) permet de simplifier la 
representation des caracteristiques. 

On ecrira en effet : 









id 


C P T 0 


1- 


[l 


7 


Cu Ahs 




IpoJ 




"" u - u 2 "" 


co 2 R 2 





(55) 



Ca _ q = qrT 0 
U p 0 AcoR p 0 AcoR 



(56) 



On obtient alors une caracteristique unique utilisable en regime subsonique (figure 17). 
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Figure 17 : Representation adimensionnelle des caracteristiques d’une turbine. 
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La puissance utile de la turbine est obtenue grace a la relation : 



P = q.Ahs.r| T 



( 57 ) 



oil r] T est le rendement global de la turbine. 

Le blocage apparait dans le distributeur et/ou la roue sous la forme d’une onde de choc 
correspondant a un nombre de Mach egal a 1’ unite. Nous examinerons ulterieurement cette 
particularite des turbines a gaz et a vapeur. 



2-2 - Methodologie de dimensionnement. 

Concernant le cahier des charges, plusieurs formes de jeux de donnees sont envisageables 
suivant le type de moteur. Nous retiendrons la forme suivante oil le nombre de donnees est 
suffisant pour resoudre le probleme : 

• La nature et les caracteristiques du fluide : 

Y' r > C = (-) (J/kg°K) (”) 

Y-l 



• Les conditionss d’entree 



• Le debit de masse 



• La vitesse utile 



• La vitesse de rotation 



po, To, p ( ) 

q 

p 



r 



N 



OU CO = 



V 



tiN 

30 



(bars), (°K), (kg/m 3 ) 
(kg/s) 

(W) 

(tr/min) 



N sera aussi elevee que possible pour satisfaire a un faible encombrement mais sera limitee 
pour des questions de tenue mecanique. 

Une premiere estimation du rendement et du nombre d’etage de la turbine est necessaire : 



r| T et n 



Calcul de la chute d’enthalpie par etage : 



et 



Ah s = — - — 
ri T qn 



( 58 ) 
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Ah = 



P 

q.n 



(59) 



Calculs concernant l’etat du fluide en sortie d’etage 

m Ahs 



2S 



c t 



■+t a 



'-'P 



P 2 =P() 

P 2 

P 2 =— =r 
r T, 



f r T' \ 
l 2S 



V T o j 



y-l 



(60) 

(61) 



(62) 

(63) 




Figure 18 : Detente dans I’etage. 

Exemple de dimensionnement global : 

On s’impose une turbine a action (annulation de la poussee axiale) et un coefficient de debit 
egal a 1'unite ou a toute autre valeur realiste (on retient souvent 0,74). On calcule 
successivement : 



r.i-°L.o 

2U 



Cu = 2U et \\t = 2 



( 64 ) 
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(65) 



U = 




R = — (taille de la machine ou rayon mo yen de l’ailetage) (67) 

co 



si la taille est trop importante, il faut augmenter co ou le nombre d’etages n. 

• cp = — — = 1 ou 0,74 
U 

d’ou Ca = cp . U (67) 

La connaissance de Ca permet de calculer l’envergure des aubages b () puis bi et b 2 (figure 19) 
dont l’importance relative par rapport au rayon R doit rester assez faible pour eviter un 
vrillage d’aubage mais egalement suffisamment elevee pour etre significativement superieure 
au jeu radial de fonctionnement. 



'blstribu t e,ur 




Figure 19 : Allure typique d’un etage a action. 

Cette situation ideale, quand elle existe, s’obtient en jouant successivement sur la vitesse de 
rotation et le nombre d’etages et eventuellement en diminuant ou augmentant cp. 



b 



o 



Q 

p 0 2n R Ca 



( 68 ) 
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avec A = 2 n R b 0 



(69) 



Triangles de vitesses : 

Les relations etablies au debut du paragraphe donnent les valeurs des angles d’entree et de 
sortie du distributeur et de la roue. 



tg cco = 0 (par hypothese) 

Cu Cu U w 

tgCC!= = = — 

1 Ca U Ca cp 

tgp, =tga 1 = 

Ca cp 

tg P, = — — = - — (par hypothese) 
Ca cp 



(71) 

(72) 

(73) 

(74) 



Definition des aubages : 

Au stade de l’avant-projet et dans le cadre du probleme inverse, un grand nombre de 

methodes sont disponibles pour definir l’ailetage du distributeur et de la roue. 

A - L’ utilisation de profils normalises correspondant a des essais systematiques en grille, tels 
que les profils NACA de la serie A3K7 pour le distributeur et de la serie B1E1I1 pour les 
ailettes de la roue (NACA TN 3802) dont les deflexions sont limites a 100° environ. 

B - La methode de l’hodographe basee sur une transformation conforme de l’ecoulement 
potentiel source-puits donnant une description correcte de l’ecoulement inter-aubages 
mais ne respectant pas la condition de periodicite. Dans cette approche, les bords 
d’attaque et de fuite ne sont pas definis avec precision. 

C - La methode des singularites permettant de definir par retouches successives la forme des 
profils sur la base d’une distribution controlee de vitesse et de pression a la surface des 
aubages. 

D - Au meme titre, les methodes numeriques dites inverses. 



Avant d’utiliser la methode des singularites ou les methodes numeriques, il est interessant 
d’ avoir recours, au prealable, a la methode A pour fixer une forme de depart (figure 15). 

Pour les etages basse-pression (BP), l'envergure de la pale devient tres importante et la 
variation radiale de la vitesse d'entrainement U modifie l'angle Pi dans des proportions 
importantes. II en resulte une evolution radiale de la section d'aubage et de son angle de 
calage (figure 20). 
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Figure 20 : Vrillcige et differentes sections d'une pale basse-pression. 
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2-3 - Blocage du debit : 



Le blocage du debit peut se produire dans l’espace interaubages du distributeur et/ou de la 
roue lorsque le nombre de Mach atteint P unite. 

L’ecoulement interaubages est determine par le champ de vitesse C/Co. On doit d’abord 
acceder aux champs de pression et de temperature pour definir le nombre de Mach local. 

Le champ de pression de la detente, supposee isentropique, est calcule par : 



Y Po 



r \ 



Y-lpc 



Po ) 



r± 

y C 2 y p 0 C ( y 

+ — = cte = — — h — — 
2 y-lpo 2 



( 80 ) 



et le champ de temperature par : 



C 2 C 

Cp T + — = cte = Cp T 0 + — y 



( 81 ) 



d’ou la vitesse du son local : 

a = Vy^T ( 82 ) 

puis le nombre de Mach local : 




Remarque : 

L’ equation (81) s’ecrit egalement en fonction de M : 



Cp T 



1 

+ 

1 1 

to 

1 


= CpT 0 


1 + ^M 0 2 


2 




2 



(83) 



y-1 2 

A P entree de l’etage, la vitesse est reduite et l’on neglige le terme — M 0 devant 1. 
On obtient alors : 
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-l+ r M 2 


T 


l P J 
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Figure 21 : Ecoulements transsoniques [2]. 




Figure 22 : Isobar es de Vecoulement trans sonique dans un distributeur[2] . 
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2-4 - Caracteristiques adimensionnelles : 

Le fonctionnement d’un etage de turbine depend de plusieurs parametres libres pouvant etre 
groupes comme suit : 

Caracteristiques du fluide : 

y, Cp, (l (viscosite dynamique) d’ou l’on peut deduire d’autres caracteristiques : 



Etat thermodynamique a l’entree : 

Pression po et temperature T 0 . Le gaz etant suppose parfait, on accede a la masse volumique 
d’entree p 0 : 

P =^~ 

Po „ 

rT o 

Geometrie de la machine : 



Elle est caracterisee essentiellement par sa taille donnee par le rayon mo yen R de l’ailetage et 
la section axiale debitante A a l’entree du distributeur (figure 19). 

Conditions de fonctionnement : 



La vitesse de rotation N (ou la vitesse angulaire CO) et la pression a la sortie pi. 

On demontre que lorsque les valeurs des 9 variables libres sont fixees (elles sont 
independantes) : 



y, C p , p, po, T 0 , R, A, co, p 2 



les variables liees suivantes sont egalement fixees au travers des caracteristiques propres de 
la machine : 



Debit massique q 
Rendement global T|t 

Ces caracteristiques sont les fonctions f et g suivantes : 



q = f (y, C p , (I, po, T 0 , R, A, (D, p 2 ) 



% = g (y, C p , p, po, T 0 , R, A, co, p 2 ) 



MACHINES A FLUIDES COMPRESSIBLES - TURBINES ET COMPRESSEURS 

Robert REY - Professeur 



28 



Le theoreme de Vashy-Buckingham indique que les fonctions f et g constitutes de 10 
variables (9 variables libres + 1 variable liee) peuvent s’ecrire sous la forme de 6 groupements 
independants et sans dimension construits a partir des 10 variables precedentes. 



Le nombre de 6 est obtenue par la difference entre le nombre de variables et le nombre de 
dimensions fondamentales intervenant dans la definition des 10 grandeurs physiques 
precedentes : dans le cas present, les dimensions fondamentales sont la longueur L, la masse 
M, le temps T et la temperature K. 

Concemant les turbines a gaz, les groupements adimensionnels sont les suivants pour la 
fonction f : 



• Le nombre de Mach debitant d’ entree ou debit reduit : 



rT A 






Y 



Po Aa o Po A 



avec a 0 = ^ yr T 0 la vitesse du son a P entree de la machine. 



• Le nombre de Mach peripherique : M LJ = — = W ^ 



l o 



V Y rT 0 



L’exposant isentropie y (deja adimensionnel) 



Le rapport geometrique de passage : 



R 

Va 



• Le rapport de detente : 




Po 

P 2 



(87) 



( 88 ) 



• Le nombre de Reynolds lie a la vitesse peripherique : Re = = (89) 

JL P r T 0 

Po 

Pour la fonction g, le parametre adimensionnel est le rendement r] T lui meme. 

Les parametres adimensionnels et les fonctions caracterisant les performances d’une turbine 
seront done les suivants : 




(90) 
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Lorsque le nombre de Reynolds est suffisamment eleve, son importance devient negligeable. 
Pour une turbine donnee (A et R constants) fonctionnant avec un fluide donne (y et r 
constants), les parametres de similitude pourront etre reduits a : 






Figure 24 : Rendement d'une turbine supersonique [2]. 



Cette representation est particulierement interessante pour les turbines supersoniques puisque 
la caracteristique debit pression se ramene a une equation simple : 




Po 



(92) 
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ANNEXE 1 



REALISATIONS TECHNOLOGIQUES 



Turboreacteurs 

Turbomoteurs 
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in - ’■ '■> 




Figure A6 : Turbomoteur utilise dans le cadre de la cogeneration. 
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Distributeur fixe Roue mobile 



Figure A7 : Sillage et couches limites du distributeur et de la roue. 




Figure A8 : Ecoulements secondaires dans le rotor. 
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Ill - LES COMPRESSEURS 



Les compresseurs se presenters sous des formes tres variees en fonction de leurs conditions 
d’ utilisation et de leurs caracteristiques de pression et de debit. 

Deux grandes classes de machines peuvent etre retenues : les compresseurs rotodynamiques 
ou turbocompresseurs axiaux et centrifuges et les compresseurs volumetriques. 

A - Les turbocompresseurs. 

Les turbocompresseurs voient un certain nombre d’application dans l’industrie mais leur 
application principale est constitute par l’aeronautique oil ils assument la compression d’air 
de la majorite des turbines a gaz (figures des Annexes 1 et 2). 

• Les compresseurs axiaux (figure 30) se distinguent par leur faible rapport de compression 
par etage et leur debit volume tres eleve. Un grand nombre de ces machines est le siege 
d’ecoulements supersoniques. 

• Les compresseurs centrifuges d’allure tres caracteristique (figure 31) ont un taux de 
compression plus eleve. 




Figure 30 : Compresseur axial (ALSTOM). 



L’empilement des etages s’effectue suivant deux principes essentiels : la construction a 
tambour et la construction a disque (figure 32). L’envergure des ailettes varie le long de 
l’ecoulement pour compenser les variations de la masse volumique du fluide et pour 
conserver a la vitesse debitante axiale une valeur constante. 
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Chambre de 
combustion £ 
flux inverse 




Prise de 
puissance 



Compresseur : 

2 Stages axiaux 
1 6tage centrifuge 



Distributee Turbine HP Turbine libre 

de turbine HP mono6tage mono<§tage 



Fteducteur 



Figure 31 : Compresseurs axiaux et centrifuge d’un turbomoteur (TURBOMECA). 



Aspiration 



Refoulement 




Piston 



(a) construction a tambour 



Aspiration 



Refoulement 




V V \ 

Disques Disque \ Piston 

fixes mobiles \ 

Labyrinthe 

(B) construction a disques 

p 0 pression d'aspiration 



Figure 32 : Construction d disques ou ci tambour [2]. 
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Les ailettes rotoriques sont rapportees sur le disque ou le tambour par divers systemes 
d’ attache (figure 33a) dont certains tres complexes en forme de sapin (figure 33b). Les 
disques statoriques sont souvent en deux parties avec plan de joint horizontal (figures de 
l’Annexe 2). 





Figure 33a : Different s systemes de fixation [2]. 




Figure 33b : Fixation des pales 
d’un compresseur axial (TURBOMECA). 



Les dispositifs d’etancheite interne et exteme conditionnent directement le rendement 
volumetrique. L’equilibrage de la poussee axiale s’effectue par un piston d’equilibrage sur 
lequel agissent de part et d’ autre les pressions de refoulement et d’ aspiration (figures 30 
et 34). 




Figure 34 : Dispositif d’equilibrage de la poussee axiale [2] 
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Pour les compresseurs terrestres, divers dispositifs mecaniques permettent de modifier 
l’orientation des aubages statoriques (figures 35) en vue de maintenir un bon rendement pour 
divers debits d’ utilisation. 




Figure 35 : Dispositif d’ orientation des aubages statoriques [2]. 
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B - Les compresseurs volumetriques. 



Les compresseurs volumetriques ont pour principe general la variation de volume d’une 
cavite entre l’entree et la sortie de la machine. Ils sont classes en deux grands types : les 
compresseurs alternatifs et les compresseurs rotatifs (figure 36). Ils se caracterisent par un 
rapport de compression pouvant atteindre des valeurs tres elevees mais sont limites en debit. 

La compression fractionnee en plusieurs etages favorise 1’ insertion de dispositifs de 
refroidissement du gaz permettant de reduire le travail moteur (figures 37). 




Figure 36 : Classification des compresseurs [2]. 




Figure 37 : Compresseur industriel a deux etages et double effet 
avec refrigeration intermediate [2 ]. 
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L’ aspect mecanique est tres soigne en particular le guidage axial, les differents dispositifs 
d’etancheite (figure 38). De nombreuses configuration de compression sont envisageables 
grace a la forme et au mode de travail des pistons (figure 39). 

Dans ces machines, les clapets d’ aspiration et de refoulement jouent un role tres important : le 
volume mort ne doit pas augmenter. Les courses des clapets doivent etre tres faibles pour 
reduire les temps d’ouverture et de fermeture mais les pertes de charge doivent etre reduites 
(figure 40). Comme pour les moteurs thermiques on peut disposer plusieurs clapets par 
cylindres. 





® vue en coupe de I'ensemble 
de compression non lubrifie 



Ch chambre de protection interdisant d la partie lubriftee de la tige de 
p6n6trer dans le cylindre 
Cl clapets enticement inoxydables 

D ddecteur d’huile arr§tant le cheminement de I'huile sur la tige du 
piston 

G garnitures PTFE d rattrapage de jeu automatique 
P piston 

R racleurs d'huile autor6glables 
SE segments d'6tanch6it6 PTFE 

SP segment porteur PTFE supportant le poids du piston 



Figure 38 : Etancheite du piston alternatif- Segments et garnitures [2]. 




Figure 39 : Differents modes de travail du 
piston [2]. 



Figure 40 : Clapet de compresseur 
alternatif [HOERBINGER] . 



MACHINES A FLUIDES COMPRESSIBLES - TURBINES ET COMPRESSEURS 

Robert REY - Professeur 



42 





Comme pour les compresseurs alternatifs, les principes utilises pour les compresseurs rotatifs 
sont innombrables (figures 41 et 42). 




Figure 41 : Compresseur ci vis [2]. 




Figure 42 : Compresseur a palettes [2]. 
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3-1 - Equations generates. 



La compression des gaz, quel que soit le type de machine, peut etre materialisee par 
P evolution des variables d’etat du fluide sur divers diagrammes. Le diagramme retenu ici est 
le diagramme T - S (figure 43). 

Pour passer du point d’entree (pression p [? temperature Ti) jusqu’a la pression de sortie p 2 , 
plusieurs evolutions sont possibles : 



• La compression isentropique (evolution adiabatique reversible ou compression calorifugee 
d’un fluide non visqueux) la temperature de sortie est egale ici a T 2 s- 

• La compression isotherme (evolution a temperature constante) necessitant un 
refroidissement tres efficace reparti tout au long de la compression. 

• La compression reelle correspondant a un compromis entre refroidissement par echange 
de chaleur avec l’exterieur et echauffement du gaz par generation d’entropie bee aux pertes 
diverses au sein de l’ecoulement. 

Pour une compression adiabatique la temperature de fin de compression T 2 sera superieure a 

T 2S ; pour une compression refroidie, T 2 sera inferieure a T 2S . 

II est interessant de comparer les travaux de compression mis en oeuvre pour chacune de ces 

evolutions. Pour cela, rappelons, P equation generale des machines de compression : 

V = f— + A ^+V (100) 

Jl p 2 




MACHINES A FLUIDES COMPRESSIBLES - TURBINES ET COMPRESSEURS 

Robert REY - Professeur 



44 



T t h : travail par unite de masse (J/kg) a foumir pour evoluer de l’etat 1 a l’etat 2. 

p : masse volumique du fluide en kg/m 3 , 

p : pression du fluide en Pa. 

C 2 

A— : variation d’energie cinetique de l’etat 1 a l’etat 2 supposee negligeable ci- 

dessous. 

if : travail de frottement du fluide (pertes de charge). 



Compression isotherme : 

En se plagant dans le cadre du fluide non visqueux (Xf = 0) et dans l’hypothese oil l’energie 
cinetique est negligeable, l’integration du terme de pression pourra s’effectuer a partir de 
l’equation d’etat des gaz parfaits : 



— = r T = rTj = cte 
P 



permettant de calculer le travail de compression : 









( 101 ) 



ou encore : 




( 102 ) 



Le travail de compression en evolution isentropique (18) permet de calculer x t h quel que soit 
1’ expos ant retenu : 




(103) 



Application numerique: 

Examinons les conditions de compression d’air depuis les conditions initiales : 

pi = 1 bar 
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Ti = 20°C 



jusqu’a la pression finale p 2 = 6 bars. 



Evolution isotherme : 



T 2 = T! = 20°C 



P, 



Pi 

rTj 



1.10 s 

287.293 



1,1 89 kg/m 3 



Evolution isentropique : 



x 0 = 1,506 J/kg 

Y= M 

7-1 






T = T 

X 2S M 



P J 



489°K = 216°C 



et x t h = 1,967 J/kg 

soit un travail superieur de 30,6% au travail isotherme. 

Cet exemple simple montre l’interet de refroidir au mieux en cours de compression. La figure 
44 montre 1’ evolution du rapport x th / x e pour diverses valeurs de k et du rapport de pression 
P 2 /P 1 . 




Figure 44 : Evolution du travail polytropique enfonction de /’ exposant k 
et du rapport de compression. 



Compression reelle : 



Dans le cas d’une compression refroidie, on fait appel a un exposant k pour modeliser la 
compression. Le travail moteur prend une forme similaire au travail isentropique oil y est 
remplace par K : 



x 



th 



Pi k 
Pi k-1 



k-l 




Pi 



(104) 
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Plus le refroidissement est efficace, plus l’exposant k est faible et plus le travail moteur est 
reduit. 



Suivant le type de refroidissement, on peut donner les valeurs moyennes suivantes pour les 
compresseurs volumetriques : 

a) - Mono - cylindre ailete : k = 1,35 

b) - Mono - cylindre ailete refroidi : k = 1,28 

par ventilateur 

c) - Bi - cylindre avec echangeur : k = 1,20 

intermediaire et ventilateur efficace 

d) - Idem avec refroidissement a l’eau : k = 1,10 

Dans les memes conditions que celles de l’exemple deja traite, on observerait les temperatures 
finales suivantes : 

Evolution isentropique : k = y = 1,4 — > T 2 = 216°C 

k = 1,28 -> T 2 = 160°C 

Evolution polytropique : < k = 1,20 — > T 2 =122°C 

k = 1,10 -> T 2 = 72°C 
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3-2 - Etage axial de compression : 



On reprendra dans ce chapitre les principales hypotheses retenues dans le cadre des etages de 
detente, a savoir : 

les forces de pesanteur sont negligeables, la compression est isentropique et le fluide se 
comporte comme un gaz parfait. 

Dans l’espace rotorique, la rothalpie reste constante. Cette propriete s’ecrit suivant plusieurs 
formes : 



W 

W 



d P C 2 

— H U Cu = cte 

P 2 

dp W 2 U 2 

— H = cte 

p 2 2 



( 105 ) 

( 106 ) 



W 2 Y p W 2 

I = Cp T + = —^~ ^ + = cte 

2 y-lp 2 



Dans l’espace statorique, il en est de meme pour l’enthalpie totale 

„ f d P C 2 
H = — + — = cte 

J p 2 



H = CpT + — = -^ P- + — = cte 

2 y-\ p 2 



( 107 ) 



( 108 ) 

( 109 ) 
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L’etage est constitue d’une roue mobile (rotor) toumant a la vitesse angulaire co et d’un 
redresseur fixe (stator) ayant pour role de ramener la vitesse absolue suivant la direction 
axiale (figure 45). 



Deux demarches bien distinctes sont a mettre en evidence : 

A - Le probleme inverse : 

II consiste a definir la geometrie d’un etage pour atteindre un certain niveau de performances : 
[rapport de pression - debit masse] pour une vitesse de rotation donnee. 

Cette demarche de "designer" necessite un certain savoir-faire et une source sure 
d’informations, correlations et resultats d’essais de toutes origines. 

Le point de depart consiste a definir a priori la distribution de travail radial et son incidence 
sur la distribution de la vitesse et de la pression (equilibre radial - voir figure 46). 

L’etablissement des champs de vitesse moyennes dans les sections 1, 2 et 3 sont effectues sur 
la base de la demarche ci-dessous ou le choix des coefficients de pression et de debit doivent 
etre effectues arbitrairement. 



Coefficient de pression : 

Ahs Cu . 






u 2 



u 



Coefficient de debit : 



Ca 
cp = — 
U 



Angles rotoriques et statoriques : 



U U 1 



w w x 

C, Ca cp 



de 0,2 a 0,6 



de 0,3 a 0,6 



tg P 2 = 



U-Cu 2 
_ Ca 



- [l-v] 

9 



tgOC 2 



Cu, 

Ca 



V 

9 



tg a 3 = 0 

Une fois ces champs de vitesse etablis, il convient de definir la forme des profils 
aerodynamiques susceptibles de les realiser. 
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r(m) 



0.100 



0.080 



0.060 




- 0.020 



0.000 



0.020 



z(m) 



0.040 



0.060 



0.080 



Figure 46 : Soufflante axiale lignes de courant en vue meridienne - Regime normal. 
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B - Le probleme direct : 



C’est l’analyse des performances d’une machine existante definie sur plans ou ayant fait 
l’objet de l’etape precedente A. L’ecoulement interaubages est d’abord determine par le 
champ de vitesse relative puis absolue C/Ci (figure 47). 

L’ equation (105) permet d’en deduire le champ de pression par la relation de St-Venant : 



Y Pi 
Y-l Pi 




IP J 



C 2 

+ - — UCu 
2 



_ Y Pi , C, 2 
Y-l p, 2 



( 110 ) 



puis le champ de temperature : 



c 2 


Ci 




Cp T + — 
2 


-UCu = Cp T. 

1 2 


(111) 


et enfin la vitesse locale du son : 


a = v/Y rT 


(112) 


et le nombre de Mach local : 


M = — 


(113) 




a 





% 



/ 




Figure 47 : Champ de vitesse absolue de V ensemble rotor-stator. 
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3-3 - Caracteristiques globales des compresseurs axiaux : 

Les caracteristiques globales se presenters generalement sous la forme du rapport de pression 
en fonction du debit masse et pour diverses vitesses de rotation. Elies sont souvent mesurees 
sur le compresseur complet c’est a dire en version multietagee (voir par exemple figure 48). 

Pour les compresseurs transsoniques, la caracteristique presente une allure verticale (figure 
49) correspondant au blocage du debit par l’apparition dans l’espace interaubages d’une onde 
de choc (figure 50). 





Figure 48 : Compresseur axial susbonique multistage du LEMFI. 
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Figure 49 : Caracteristique globule pression-debit. 




Figure 50 : Blocage du debit dans Vespace interaubages rotorique [ONERA]. 
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Comme pour les turbines, on demontre que les parametres adimensionnels se ramenent a : 




( 114 ) 



La ligne de pompage limite 1’ utilisation du compresseur en debit partiel : fonctionnement 
instable sur le plan vibratoire lie a des decollements toumants a l’entree du rotor. 
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3-4 - Compresseurs volumetriques alternatifs : 



Quelques elements importants seront abordes ici, ils concernent tous les diverses pertes des 
compresseurs volumetriques. 

A - Coefficient de remplissage et espace mort : 

Observons le fonctionnement du compresseur sur la figure 51 suivante ou les clapets sont 
supposes parfaits. Les temps d’ouverture et de fermeture des clapets conditionnent 
directement le rendement volumetrique du compresseur : en phase de refoulement, une partie 
du debit aspire est refoule dans 1’ aspiration (retard de fermeture du clapet d’ aspiration). 

En phase d’aspiration, une partie du gaz deja refoule est reaspire au travers du clapet de 
refoulement (retard de fermeture du clapet de refoulement). 




Figure 51 : Cycle de compression enfonction du balayage du piston. 



On definit sur le diagramme precedent les differentes phases : 

BC : compression 

CD : debit a pression constante 
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DA : detente du gaz contenu dans l’espace mort 

AB : aspiration a pression constante 

On definit egalement les volumes suivants : 

Vth : volume balaye par les pistons 

Va : volume reellement aspire 



Le coefficient de remplissage est defini par le rapport : 



Bv = 



Va 

Vth 



On demontre que celui-ci prend la forme : 



1 1 / k 






1-g 



-1 



oil g represente l’espace mort. 



( 115 ) 



(116) 



Pour — = 6.32 et k = 1.3 on observe les valeurs suivantes : 



Pi 

G 

G 

G 

G 

G 



= 1% 


r| v = 96,8 


= 2% 


r| v = 93,7 


= 3% 


r| v = 90,6 


= 4% 


r| v = 87,5 


= 5% 


T| v = 84,3 



Pour un meme compresseur, le passage de 1% a 5% d’espace mort entraine une baisse de 
debit de l’ordre de 15%. 



Pression de refoulement maximale : 



C’est la valeur pour laquelle le debit du compresseur s’annule. On obtient cette condition 
limite a partir de la relation precedente lorsque le coefficient de remplissage s’annule : 



soit : 



0 = 1-g 





p p 


1/k 




P 2 max 


— 1 




Pi J 





P 2max 




Pour pi = 1 bar et k = 1,3 on obtient : 



= 1% 


P 2 max = 403 bars 


= 2% 


P2max= 166 bars 


= 3% 


P2 max = 99 bars 


= 4% 


P2 max = 69 bars 


= 5% 


P 2 max =52 bars 
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B - Influence des pertes de charges : 



Dans un compresseur volumetrique les pertes de charges sont de plusieurs nature mais 
essentiellement liees aux changement rapide de section et de direction (clapets) ainsi qu’aux 
frottement dans les tuyauteries de liaison et eventuellement l’echangeur thermique. Ces 
dernieres pertes sont directement liees aux vitesses d’ecoulement et done aux sections de 
passages offertes. Si l’on suppose que les pertes dans les clapets d’aspiration et de 
refoulement sont egales a Ap, le travail moteur prendra la forme suivante : 



T _Pi k 
Pi k-1 



P 2 + A P 
Pi - Ap 




(117) 



Toujours pour le meme exemple numerique avec k = 1,35 et Ap = 0,3 bar, on obtient les 
valeurs suivantes : 

sans pertes : x t h=l,98.10 5 J/kg 

avec pertes : x = 2,56.10 5 J/kg 



soit une augmentation de 29,3% du travail a foumir au fluide. 



C - Influence du refroidissement : 



Le compresseur est constitue de deux etages en serie. Le refroidissement est assure par des 
ailettes equipant les corps de cylindre ainsi que par un refroidissement intermediate a 
serpentin situe entre la sortie du premier etage et L entree du second. L’ensemble est refroidi a 
Pair par un ventilateur a roue ouverte (figure 52). 




Un dispositif experimental est mis en place et permet de relever les valeurs des differentes 
pressions et temperatures. 

Ces essais ont pour but de definir les exposants polytropiques k de chacun des deux etages 
ainsi que de V ensemble du compresseur. 
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Les resultats sont donnes dans le tableau numerique ci-dessous pour diverses valeurs de la 
pression de refoulement : 



Po 

bars 


To 

°C 


Pi 

bars 


Ti 

°C 


P2 

bars 


t 2 

°c 


P 3 

bars 


t 3 

°c 


1 


20 


4,8 


166,3 


4,25 


54 


6,7 


75 


1 


20 


5 


174,4 


4,40 


56 


9,8 


101 


1 


20 


5,5 


184 


5,00 


56,8 


20,8 


146 


1 


20 


6,0 


190 


5,40 


56,3 


31,5 


168 



La relation general polytropique (21) appliquee au l er etage, au second etage et a l’ensemble 
du compresseur donne les valeurs suivantes : 

l er etage : kj = 1,35 

2® me etage : k 2 =l,19 

ensemble : k = 1,12 

Dans l’absolu, ces valeurs sont acceptables et traduisent un refroidissement glogalement tres 
convenable hormis cependant le premier etage. On notera l’excellente du refroidissement 
intermediate. 



D - Determination des pertes de clapets : 

Les pertes de clapet sont calculees a partir des resultats d’essais effectues a vide et obtenus en 
debranchant toutes les tuyauteries entrant et sortant des 2 cylindres puis en mesurant les 
temperatures correspondantes. On trouve : 

l er etage : temperature de sortie = 60°C 

2 eme etage : temperature de sortie = 3 8°C 

grace a ces valeurs on peut definir Ap du l er et 2 eme etage par la relation : 



T. 

T, 




k-l 

k 



ou k prend les valeurs issues du paragraphe precedent. 
On a ici : 



p 2 = po + Ap 

Pi = Po - Ap 
T, = 20°C = 293°K 
po = 1 bar 

On obtient pour le l er cylindre Api = 0,242 bar et pour le 2 eme cylindre Ap 2 = 0,185 bar. 
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ANNEXE 2 



REALISATIONS TECHNOLOGIQUES 
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Figure B1 : Compresseurs en cours de montage (SULZER). 
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Figure B2 : Compresseurs centrifuges et axial (RATEAU). 
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Figure B3 : Compresseurs centrifuge el axial - Dessin 3D. 
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R.l. - Produit vectoriel 



— ^ ^ 

Soient deux vecteurs V| et V 2 definis par leurs composantes respectives (coordonnees 
cartesiennes) : 





U 1 




u 2 


Vt = 


Vl 


et V 2 = 


v 2 




Wj 




w 2 



On definit le produit vectoriel par : 



V-AvU 



V 1 W 2 - WjV 2 
W}U 2 - U|W 2 

UlV 2 - Vjllj 



— > — > 

Si Vj et V 2 sont colineaires. 



Vj A V 2 = O 



R.2. - Produit scalaire 



— > — > 



Le produit scalaire de V, par V 2 est defini par 



— > — > 



V|. V 2 = UjU 2 + + WjW 2 = VjV 2 cosa 



— > — > 



— > — > 



Si V, et V 2 sont perpendiculaires, . V 2 = O 



R.3. - Gradient d'une fonction scalaire cp 

La fonction (p (x, y, z) a pour gradient le vecteur suivant : 
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grad cp = 



dcp 





-4 

grad cp = 



dtp 

dr 

d£ 

rd0 

dcp 

dz 



coordonnees cartesiennes 



coordonnees cylindriques 



R.4. - Divergence 



— ^ 

Soit le vecteur V 



u 

v , on definit sa divergence par le scalaire 



w 



div V 



du dv dw 
dx dy dz 



(coordonnees cartesiennes) 



R.5. - Rotationnel 



— ^ 

Le vecteur V a pour rotationnel le vecteur : 



dw 


dv 


dy 


dz 


du 


dw 


dz 


dx 


dv 


du 


dx 


dy 



R.6. - Laplacien 

Laplacien d'une fonction cp(x, y, z) : 

d 2 cp d 2 cp d 2 cp 
V2<P = dx* + dy* + dz* 



et en coordonnees cylindriques : 

V2fft _ J_ d^cp d^cp 

9 dr 2 r dr r 2 d0 2 dz 2 



MACHINES A FLUIDES COMPRESSIBLES - TURBINES ET COMPRESSEURS 

Robert REY - Professeur 



66 




R.7. - Theoreme de Green-Ostrogradski 



|j s n . AdS = Jj| v div A . dV 



R.8. - Differentielle totale 

Soit (|) une fonction des coordonnees et du temps : 

4> = 4>(x,y,z,t) 



La differentielle totale de 4> sera donnee par : 

,, 9ct) 3d) 3d) 3d) 

d(|) = — dx + — dy + — dz + — dt 

3x 3y 3z 3t 



Inversement, la fonction K = K , dx + K 2 dy + K 3 dz + K 4 dt est une differentielle totale 
exacte egale a dU si l'on a simultanement : 

3U au 3U 3U 

JA.1 , Pvp — , J\o — , IV/l — 

1 dx - 3y 3 3z 4 3t 

R.14. - Mouvement relatif 



Ecrivons l'irrotationnalite de l'ecoulement absolu : 



3 

rot C = 0 



En coordonnees cylindriques, cette relation vectorielle s'ecrit : 



1 3Ca 3Cu 
r 30 3z 



3 

rot C = 



3Cr 3Ca 
3z 3r 



1 3(rCu) 1 3Cr 
r 3r r 30 



(1) 
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En remplagant chaque composante de la vitesse absolue par se valeur en fonction des 
composantes relatives et d'entrainement : 





0 


Wr 


c = tJ + w = 


cor + 


Wu = 




0 


Wa 



Cr = Wr 
Cu = cor + Wu 
Ca = Wa 



Le systeme (1) s'ecrit : 



> 

rot C 



1 5Wa 


5Wu 


r 5r 


5z 


5Wr 


5Wa 


5z 


5r 


1 5(rWu) 

Oz-.-v 1 


r 

r 


5r 



= 0 



1 5Wr 
r 50 



= 0 



En observant que le rotationnel de U s'ecrit : 



> 

rot tJ 



0 

0 

2co 



Le systeme (2) peut s'ecrire : 



(2) 



> > 

rot C - rot W + 2 co 



ou encore : 



> 

rot W 



1 5Wa 5Wu 
r 5r 5z 

r)Wr _ f)Wa 
dz r)r 

1 5(rWu) 1 5Wr 
r 5r r 50 ^ 



Le mouvement relatif est done rotationnel. 



(3) 
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